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Abstract 

This work examines the magnetic order and spin dynamics of a double-exchange model with com- 
peting ferromagnetic and antiferromagnetic Heisenberg interactions between the local moments. 
The Heisenberg interactions are periodically arranged in a Villain configuration in two dimensions 
with nearest-neighbor, ferromagnetic coupling J and antiferromagnetic coupling —nJ. This model 
is solved at zero temperature by performing a 1/yS expansion in the rotated reference frame of 
each local moment. When n exceeds a critical value, the ground state is a magnetically frustrated, 
canted antiferromagnet. With increasing hopping energy t or magnetic field B, the local moments 
become aligned and the ferromagnetic phase is stabilized above critical values of t or B. In the 
canted phase, a charge-density wave forms because the electrons prefer to sit on lines of sites that 
are coupled ferromagnetically. Due to a change in the topology of the Fermi surface from closed 
to open, phase separation occurs in a narrow range of parameters in the canted phase. In zero 
field, the long-wavelength spin waves are isotropic in the region of phase separation. Whereas the 
average spin-wave stiffness in the canted phase increases with t or 77, it exhibits a more complicated 
dependence on field. This work strongly suggests that the jump in the spin-wave stiffness observed 
in Pri_ x Caa;Mn03 with 0.3 < x < 0.4 at a field of 3 T is caused by the derealization of the 
electrons rather than by the alignment of the antiferromagnetic regions. 



I. INTRODUCTION 



The persistence of antiferromagnetic (AFM) short-range order in the ferromagnetic (FM) 
phase of the manganites has been recognized for many years Q]. In metallic manganites 
like Lao.7Cao.3Mn03 that contain a preponderance of AFM-coupled polaronic regions 
the Curie temperature Tq is suppressed but the magnetoresistance is strongly enhanced. 
Close to but below Tq, the spin dynamics of Lao.yCao.sMnOa contains both a propagating 
spin-wave (SW) branch from the FM regions and a diffusive component from polaronic 
regions with suppressed FM interactions 4|. The low-temperature insulating phase of the 



manganite Prx^Ca^MnOs with 0.3 < x < 0.4 was originally believed p,lal2| to be a canted 
AFM (CAF) but may actually contain regions with both FM and AFM short-range order 
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131 ] - When an applied field B exceeds about 3 T, the resistivity drops by 



several orders of magnitude p, the AFM regions shrink |12j, and the SW stiffness D sw jumps 
by a factor of 3 [11]. Despite the central role played by AFM interactions in the manganites, 
little is known theoretically about how they affect the propagating SW dynamics. As first 
shown by Anderson and Hasegawa Q|, the effective electron-hopping energy between two 
local moments making a relative angle G is proportional to cos 0/2 in the limit of large 
Hund's coupling. So AFM interactions may sup press the contribution of electron-mediated 
double-exchange (DE) to the SW dynamics [l5|. By aligning the local moments, a magnetic 



field or electron hopping will alter the DE contribution to the SW dynamics. This paper 
examines the effects of AFM interactions on the ground-state properties and SW dynamics 
of eiectrons coupled to the ,ca, _ of a Villain tnodel UUW- 

The generalized Villain model is one of the simplest periodic models to exhibit magnetic 
frustration. As described in Fig. 1(a), the local moments Sj are coupled by the FM interaction 
J along the y direction and by either the FM interaction J or the AFM interaction —rjj 
along the x direction. The CAF phase is stable when rj exceeds the critical value r] c , which 
is 1/3 when B = Bz = but increases as B increases. Due to the different environments of 
the a and b sites, the angle 8b at the b sites is always larger than 8 a at the a sites, as shown 
in Fig. 1(b). 

In the hybrid model considered here, the Heisenberg interactions between the local mo- 
ments are in the Villain configuration while electrons with density p = 1 — x are FM coupled 
to the local moments by Hund's coupling Jh and hop between neighboring sites with energy 
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FIG. 1: (a) The generalized Villain model with Heisenberg couplings J or — i]J, (b) the local 
moments in the xz plane subtend angles 9 a and 9b with the z axis, and (c) the electron spins also 
lie in the xz plane but subtend angles ip a < 9 a and tpb < 9b with the z axis. 

t. The DEV model (so called because it combines the DE and Villain models) provides 
several advantag basis for understanding the effects of AFM interactions and non- 

collinearity on the spin dynamics. First, except in a narrow range of parameters, a homo- 
geneous CAF phase is stable against phase separation when the AFM control parameter rj 
exceeds r] c . By contrast, the well-studied hybrid model with AFM interactions — J between 
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mome nts p hase separates before the AFM exchange J is large enough 



221 23| . This phase instability is caused by a "site-local continuous 



l24j that is absent in the DEV model. Second, unlike the case in a 
hybrid model with AFM exchange only, the ground state of the DEV model contains a FM 
component even when t = and B = 0. So it can be used to evaluate the change in SW 
stiffness D sw as the electrons become mobile. Third, because the DEV model contains both 
FM and AFM Heisenberg interactions, it can be used to study insulating manganites like 
P r o.66Cao.34Mn03, where the AFM interactions are produced by superexchange and the FM 
interactions by short-ranged orbital order j3, 13 . 
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For simplicity, we have constructed a model that is translationally symmetric in two 
dimensions. Since this model is solved at zero temperature, the qualitative results will be 
unchanged in three dimensions. More problematically, the AFM interactions are arranged 
periodically rather than in clusters. In the low-temperature phase of Pr .67Cao.33Mn03, the 
FM interactions may be confined to two-dimensional sheets with widths of roughly 25A 



in a "red cabbage" structure |l2j, |13|. Neutron-scattering results on the other hand, 
suggest that the FM clusters in the insulating phase are about 40A is diameter. So for 
wavelengths much longer than 40A, the SW's will average over the FM and AFM regions. 
Hence, the DEV model will provide qualitatively accurate predictions for the average SW 
stiffness = (D* w + Z^f w )/2 defined in the long-wavelength limit. 
The Hamiltonian of the DEV model is given by 

H = -tJ2 + c] a c ia ) - 2 J H £ s, • S, - £ J v St ■S J -B2 S IZJ (1) 

where c\ a and c ia are the creation and destruction operators for an electron with spin a 
at site i, Sj = (l/2)c\ a a a pc i Q is the electronic spin, and Sj is the spin of the local moment 
with magnitude S. Nearest-neighbor Heisenberg interactions Jij take the values J (FM 
interaction) or —i]J (AFM interaction), as described in Fig. 1(a). This model is solved at 
zero temperature by expanding the Hamiltonian in powers of 1/ y/S. To guarantee that the 
contributions to the SW frequencies from hopping and from the Heisenberg interactions are 
of the same order in 1/y/S, t is considered to be of the same order in l/y/S as JhS, JS 2 
and BS. Hence, the dimensionless parameters of our model are t' = t/JS 2 , rj, B' = B/JS, 
and Jh/ JS. To lowest order in 1/S, the magnetic field B only couples to the local moments 
and not to the electrons. While the theory developed below can be extended to treat all 
values of the Hund's coupling, we shall for simplicity consider the limit of large J^S or in 
dimensionless terms, Ju/JS ^> 1 and J#S/t 3> 1. 

This paper is divided into five sections. The Villain model is discussed in some detail in 
Section II, where we provide new results for the SW stiffness. The ground-state properties of 
the DEV model are presented in Section III. In Section IV, we evaluate the SW frequencies 
of the DEV model. Section V contains a discussion and summary. A short version of this 



work 



is m press. 
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II. GENERALIZED VILLAIN MODEL 



This section presents, for the first time, a Holstein-Primakoff expansion for the generalized 
Villain model. In an equivalent approach, Saslow and Erwin fljj| numerically evaluated the 
mode frequencies by linearizing the equations of motion for the spins. However, a formal 
Holstein-Primakoff expansion is required to lay the foundation for the solution of the full 
DEV model in Section IV. 

The Hamiltonian of the generalized Villain model is given by Eq.Q with t = and 
Jh = 0. The spin dynamics is immensely simplified in the rotated reference frame for 
each spin: S, = U^ ot Si, where U^ ot is the unitary rotation matrix for site i. A Holstein- 
Primakoff expansion is performed within each rotated reference frame: §i z = S — ajaj, 
Six + iSiy = V^Scii, and S ix — iSi y = y/2Sa\. Then the zeroth-order term (in powers of 
1 / vS) in H v can be written as 

E h = ^NJS 2 1 - cos 29 a + rj cos 29 b - 2 cos(9 a - 9 b ) - B' (cos 9 a + cos 6 ) j , (2) 

which is of order JS 2 . 

Minimizing with respect to 9 a and yields the relations 

sin 29 a + sin(# a - 9 b ) + -B' sin 9 a = 0, (3) 

- r] sin 29 b - sm(9 a - 9 b ) + ~B' sin 9 b = 0. (4) 

In zero field, it is easy to show that 9 b = 39 a for all 77. The equilibrium angles are plotted 
versus r] for several different values of B' ranging from to 3 in the inset to Fig. 2. In the 
limit of large rj with B' = 0, 9 b —>■ n/2 and 9 a — > n/6. For nonzero field, 9 b still approaches 
7r/2 but 9 a approaches an angle smaller than n/6. In Fig. 3(a), we plot the equilibrium angles 
versus B' for several values of rj. 

After linearizing Eqs.(|3*]l and (@J), we find that the phase boundary between the CAF and 
FM phases satisfies the relation 



B' -2 V + 4-2^/(7] + 1) 2 + 1 = 0, (5) 

which was first obtained by Gabay et al. While r] c = 1/3 in zero field, r] c increases with 
the field B' as seen in the inset to Fig.2 and in Fig. 3(a). 
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FIG. 2: The SW stiffnesses versus rj for t' = and various values of the field B' . In the inset are 
plotted the equilibrium angles 6 a and 9 b versus rj for the same set of fields. 

Expanded as H v = Eh + H v \ + H v2 + . . . in powers of 1/ y/S, the first-order term H v \ 
vanishes provided that the angles 9 a and 9 b satisfy Eqs. (J2J) and (@J). In terms of the Fourier- 
transformed spin operators and a£ ^ on the r = a or b sublattice, the second-order term 



H v2 can be written as 

rj , c pf (r)t (s) . / (r) (s) . (r)f (s)A D (r,s) \ fa s. 

H v2 = JS 2^j«k «k + ( a -k«k + a -k a k J^k h ( 6 ) 

A|f a) = 2 cos 2fl a + 2 cos(fl a - 6 ) - 2 cos 2 fl a cos ^ + B' cos a , (7) 

A (a,b) = A( b,a) = _ 2 cos 2((0 q _ 6b)/2) CQS ^ (8) 

A% ' b) = 2 cos(# a - 6 b ) - 2r] cos 26 h + 2?y cos 2 6 b cos ^ + 5' cos 6 b , (9) 

5jf' a) = sin 2 fl a cos£; x , (10) 

B (a,b) = ^(M = sin 2 (( ^ _ 9b)/2) CQsk ^ (n) 

Bg' b) = -r}sm 2 d b cosk x , (12) 
where the lattice constant is set to one. 
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Ri = 4{4 a ' a)2 + A {b > b)2 - A(Bt' a)2 + Bi b > b)2 )}(A^ b)2 - B^ h 

2 



+ {At a)2 - A {b ' b)2 - 4(4 a ' a)2 - B {b ' b)2 )} + 8(4°' a) 4 6 ' 6) + AB^ a) B { £' b) ) (14) 
x(4 a ' fc)2 + B { ^ b)2 ) - 324 a ' 6) 4 a ' 6) (4 a ' a) 4 fc ' 6) + A (b > b) B^ a) ). 

The SW frequency to^ = satisfies the condition cUk=o = B. Our results for agree 
with the mode frequencies numerically evaluated by Saslow and Erwin Q]. In the long- 
wavelength limit, the SW stiffnesses in the x and y directions are obtained from the expres- 
sion lim k ^ uj k = B + ir K /,•; + D^k 2 y . 

In the FM phase, the SW frequency can be solved analytically: 



cj k = B + JS{3 -77 + (77 - 1) cos£; x .) - JSyJ(l + r/) 2 (l - cos k x ) 2 + 4 cos 2 k y . (15) 

So for k x = 0, = B + 2JS(1 — cosk y ) is independent of rj. The SW stiffnesses in the FM 
phase are given by the simple results D* w = (JS/2)(1 — rj) and D^ w = JS, independent of 
field. The FM phase becomes unstable when cuq = 0, where Q = (jr, 0) is the AFM Bragg 
vector. This yields the same condition for the CAF-FM phase boundary as Eq.©. 

In the CAF phase, analytic results for the SW stiffnesses were found only when B = 0: 



D: w = JSV2 V , 



\ V^ + 1 ' 



(16) 



V-V^t (17) 



Dl w = JSV2, 

When rj — > 00, D* w — > JSy/2rj and _D^ W — > JS^2/r/. For B > 0, -D*f w tends to a nonzero 
limit but D* w still diverges as rj — > 00. So the average SW stiffness always diverges when 
77 — > 00 regardless of the field. 

As shown in Fig. 2, Z};f w is a linearly decreasing function of rj in the FM phase below 
r] c . In the CAF phase above rj c , increases and decreases with rj. Long-wavelength 
SWs become isotropic when D* w and D*! w cross. For B = 0, = D^ w when rj — 1. This 
crossing point moves to progressively larger values of 77 with increasing field. 

The SW stiffnesses are plotted versus field in Fig. 3(b) for four different values of 77. In 
the FM phase above B c , the stiffnesses are independent of field. But in the CAF phase 
below B c , the dependence on field is more complex. For 77 > 1, D* w > D% w at zero field and 
the stiffnesses cross as the field increases. If the SW stiffness in the a direction is defined 
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(k x /7C,k y /7C) 

FIG. 4: The SW frequencies for p = 0.66, rj = 1, t' = 0, and various fields S'. 

at fixed k a by D^ w (k a ,B) = (o^ — B)/k^, then D° w (k a ,B) increases rapidly near the field 
B* = D*„(k n , 0)k^ as shown in the inset to Fig. 3(b) for fc Q = 0.0157T. As discussed in detail 
elsewhere |29j, this behavior is typical of any CAF with a quadratic SW dispersion. Since 
the drop in D* w is then steeper than the rise in D^ w , D™ decreases with field for B 3> B*. 

Surprisingly, Fig. 3(b) indicates that the stiffness in the x direction becomes negative for 
sufficiently large rj. To understand this behavior, we have plotted the SW frequencies versus 
k for r\ = 1 and several fields in Fig.4. Since uq = and u = B in the CAF phase, the 
SW stiffness in the x direction must decrease as the field B becomes comparable to JS. 
When 7] = 1, D x svl = for all fields above B c = 2(^5 - 1) JS w 2.472JS. For rj > 1 and 
B > B c , D* w < in the FM phase. Due to the reduced symmetry of the CAF phase, the 
first Brillouin zone extends from —n to n along k x but from —rr/2 to n/2 along k y j^. 

The SW frequencies are plotted versus k for B = and various values of rj in Fig. 5. 
When 7] = —1, all of the Heisenberg interactions equal J and the SWs are isotropic. For 
rj < r] c = 1/3 in the FM phase, the SW frequencies are independent of rj along k y but not 
along k x , as implied by Eq. (|T5|) . With increasing rj in the CAF phase, the SW stiffness 
increases along the x direction but decreases along the y direction, as predicted by Eqs. fjlfij) 
and (|17jl and shown in Fig.2. Also notice that the SW velocity at Q is an increasing function 
of rj in the CAF phase. 

Both Figs. 4 and 5 indicate that the SW velocity at Q softens as the CAF phase becomes 
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FIG. 5: The SW frequencies for B' = 0, t' = 0, and several values of ij. 

unstable. In the FM phase, dw^/ dk x vanishes at Q for all B > B c or rj < r\ c . The softening 
of the SW velocity at the CAF-FM phase boundary is a quite general result discussed by 
Roman and Soto j3l|. 

III. GROUND-STATE PROPERTIES OF THE DEV MODEL 

In this section, we discuss the ground-state properties of the DEV model, which is con- 
structed by placing electrons with density p on the Villain lattice. An electron on site % is 
coupled to the local moment at that site by Hund's coupling — 2JnSj • Sj and is allowed to 
hop to neighboring sites with hopping energy t, as specified by the DEV Hamiltonian of 
Eq.flU). 

With the Fermion creation and destruction operators and c£2 defined in the reference 
frames of the local moments, the band Hamiltonian of the electrons can be written 

H b = J2 £ + KW)^c% (is) 

k i,j a,/3 

where Cy/ a = {cj^, c^.n a , cj^, c^+q a } defines the ij subspace. The sum over k is restricted 
to the first Brillouin zone and a = ±1 corresponds to spin up or down in the local reference 
frames. The zeroth-order band Hamiltonian includes just the Hund's coupling: H^ a0 = 
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—JuS5ij<Ta8- The potential V%{\&) is smaller by t/JnS and has matrix elements 



Vg(k) = -V a 2 |(k) = -2t cos k x cos 9 a 5 a/3 , (19) 

F a 3 |(k) = -O k ) = ^tcos^cos^, (20) 

Kg(k) = -Vg(k) = -2tf cos sin a ^, (21) 

7$(k) = -Kg(k) = -2ti cos k x sin (22) 

^|(k) = O k ) = O k ) = V$<V = -2tcosk y cos((6 a - 9 b )/2)5 aP , (23) 

V^(k) = V$(k) = -KS(k) = -VS(k) = 2tzcos^sin((^ - 9 b )/2)a%, (24) 



Notice that the a v af3 terms couple the up and down spin states. Like E h of Eq.(J2J), H b is also 
of order JS 2 . 

The potential V^g(k) is treated within degenerate perturbation theory. Second-order 
perturbation theory with corrections of order t/JnS will be required to obtain the SW 
frequencies in the next section. But to order (t/JnS) or to first order in the potential, the 
spin up and down subspaces decouple and the band Hamiltonian H b is easily transformed 
into the diagonal form H b = T,k,a,r e ka^ka 42 by the rotations cjfj = «k <4c! + u^d^l and 



c ka - % «k Q _ % «k Q wnere e ka - -J H ba + e k , 

u (a)2 _ ^ _ ^(6)2 = i/l + — (cOS# a - COS 6 ) COs/c x |, (25) 
^ L ^k ' 

= — £ cos A^. (cos a + cos 9 b) =F ^k, (26) 
u>k = y (cos 9 a — cos 9 b \ cos 2 k x + 4cos 2 ((6> a — b )/2) cos 2 (27) 



The upper and lower signs in e k refer to the r = a and b bands, respectively. With H b in 
diagonal form, E b = (H b ) is easy to evaluate. 

It is straightforward to minimize the zeroth-order energy E = + E b with respect to 
the angles 9 a and 6 b in the limit of large J^S. The relations that generalize Eqs.@5J) and (jU 

are 

sin2£ a + sin(£ a - 9 b ) + -B'sin9 a + Ej^J-/^) + j = 0, (28) 

- ,-n2ft - sin( S „ - « + - B <sin* t + _£{Jk. /( £>) + ^M" )} = 0, (29) 
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where /(e) = 0(/2 — e) is the Fermi function at T = and p, = ji — J^S sgn(p — 1) is 
the shifted chemical potential. For a fixed rj and B', the equilibrium angles decrease with 
increasing t'. The phase boundary between the CAF and FM phases is now given by the 
condition 

B' - 2 V + 4 + 3K/AJS 2 - 2^(1 + rj) 2 + 1 + K/AJS 2 + {K/8JS 2 ) 2 = 0, (30) 

where K = -((e£°) + (ef))/ 2 > is the average kinetic energy of the electrons in the 
FM phase. This reduces to Eq.(|5|) when K = 0. For p = 0.66, r\ = 2, and B' = 0, 
the dependence of the equilibrium angles on if is plotted in Fig.6(a). Also shown is the 
average spin M = S(cos8 a + cos^)/2 of the local moments. When if = 10 and B' = 0, the 
equilibrium angles and M/S are plotted versus rj in Fig. 7(a). Their dependence on field B' 
is plotted in Fig. 8(a) for if — 3 and rj = 2. 

Surprisingly, the electronic occupation of the a and b sites on the Villain lattice are 
different. For p < 1, the occupancies of the a and 6 sublattices are given by 

na = ^EH a)2 M a) )+^ )2 f^% (31) 
iV k 

n,= ^EK ) V(4 a) )+«L 6)2 M a) )}, (32) 

From the relation m^ 2 + lij^ 2 = 1, it follows that (n a + n&)/2 = p is just the average number 
of electrons per site. 

In the CAF phase, electrons prefer to sit on the a sites of the Villain lattice. The fraction 
fa — n a/2p > 1/2 of such electrons is plotted in Figs. (6-8). For rj — 2 and B = in Fig. 6 (a), 
f a has a maximum of 0.574 as t' — > and approaches 1/2 as if — *■ 1f c ~ 13.2. Similar behavior 
is found in Figs. 7(a) and 8(a), where f a is shown to be an increasing function of rj and a 
decreasing function of B' . This behavior is easy to understand: the largest angles between 
neighboring spins are along the x axis between b sites with angles differing by 29 . When an 
electron hops onto a b site, it cannot easily hop to other b sites and so quickly moves onto a 
neighboring a site, where it can readily travel between other a sites with angular difference 
29 a -C 29b- Hence, the non-collinearity of the local moments quite naturally produces a 
charge- density wave (CDW) with a substantial amplitude. As rj — >• oo, 9b approaches n/2 
and f a approaches 0.593. So even when the electrons are unable to hop between sites on 
the b sublattice (since the angles on neighboring b sites differ by n), roughly 40% of the 
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FIG. 6: (a) The angles 6 a and 9b of the local moments, the total local magnetization M/S, the 
fraction f a of the electrons on the a sites, and (b) the spin-wave stiffnesses (both in the x and y 
directions and their average) for p = 0.66, rj = 2, and B' = versus i'. 



electrons can still be found on b sites at any one time. Because neither the CAF nor FM 
densities-of-states contain a gap, both phases are metallic within the DEV model for t' > 0. 
Due to short-range orbital ordering, a CDW with the same period as the one predicted 
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FIG. 7: (a) The angles 6 a and 9b of the local moments, the total local magnetization M/S, the 
fraction f a of the electrons on the a sites, and (b) the spin-wave stiffnesses for p = 0.66, if = 10, 
and B' = versus rj. 

here has in fact been observed in the AFM regions of Pr .7Cao.3Mn03 [7 , 32 1 . However, the 
observed charge ordering is perfect: all of the Mn 3+ ions lie on one sublattice and all of the 
Mn 4+ ions lie on the other. Such perfect charge ordering is never achieved within the DEV 
model. 
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FIG. 8: (a) The angles 6 a and Ob of the local moments, the total local magnetization M/S, the 
fraction f a of the electrons on the a sites, and (b) the spin-wave stiffnesses for p = 0.66, t' = 3, 
and rj = 2 versus B' . 



Another surprising result is that phase separation occurs within a narrow range of t'. 
In a plot of filling p versus chemical potential //, phase separation appears as a jump Ap 
in p(fj). For the parameters rj — 3 and B = in Fig. 9, Ap reaches a maximum of about 
0.0028 when t' rts 10.0 and shrinks as t' increases. If p is fixed at 0.66, then phase separation 
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2 

|i/JS 

FIG. 9: The filling p versus chemical potential jl for various values of t' , rj = 3, and B' = 0. Phase 
separation appears as discontinuities in p. 

occurs within a very narrow range of t' between about 9.98 and 10.02. For fixed p, phase 
separation appears as jumps in the equilibrium angles 9 r and in the electron fraction f a , as 
seen in Figs. 6 (a) and 7(a) for zero field. 

Like the Pomeranchuk instability 0, [mJ in the two-dimensional Hubbard model, the 
phase instability in the DEV model occurs close to but slightly above the Van Hove filling 
and is marked by a change in Fermi surface (FS) topology from closed to open. In the 
usual Pomeranchuk instability, however, the change in FS topology spontaneously breaks 
the square symmetry of the lattice. Square symmmetry is already broken in the DEV model 
by the Villain arrangement of the Heisenberg interactions. The FS of the DEV model is 
sketched in Fig. 10 for rj — 3, B — 0, p — 0.66, and for values of t' on either side of the 
phase-separated range. For t' just above 10.0, the extra electrons in the neck of the a FS 
around k = are offset by the holes in the h FS around k = (7r,7r/2). Although phase 
separation occurs for any t' around some value of the filling, it becomes significant only if 
the b FS is already present when the necks in the a FS develop. For the parameters in Fig.9, 
this requires that t' > 9.9. 
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FIG. 10: The FS is plotted for rj = 3, B' = 0, p = 0.66, and t' = 9.8, 10.2, 12, 16, or 21.2. For 
if > 10, the a FS develops necks around k = 0. 

A magnetic field very quickly narrows and then eliminates the region of phase separation. 
For t' = 3 and r\ — 2, phase separation does not occur with increasing field. The small kinks 
in the equilibrium angles and electron fraction seen in Fig.8(a) correspond to points of Van 
Hove filling (where the necks in the a FS first appear) but not to a phase-separated region. 



IV. SPIN DYNAMICS OF THE DEV MODEL 



In this section, we evaluate the SW frequencies for the DEV model by following the 
general formalism developed by Golosov 23]. While introducing some new physics, the 
presence of non-collinear spins also complicates things a bit. 

When JuS/t is finite, the relationships given above for o£i 111 terms of cf£a must contain 
admixtures of the opposite-spin terms (t/ JnS)d^_Q _ a . As sketched in Fig. 1(c), this implies 
that the equilibrium angles ipr of the electrons are smaller than the angles 9 r of the local 
moments because the electrons try to align their spins as much as possible. For large 
Hund's coupling, T — ip r (xt/ J^S and the electrons always exert a small torque on the local 
moments. In the local reference frame of site i on the r sublattice, (si X ) = (n r /2) sm(6 r — ip r ) 
so these new terms produce a correction to the Hund's coupling — JhSj ■ s, that survives in 



17 



the JuS — > oo limit. Indeed, these torque terms are required to obtain sensible results for 
the SW frequencies. 

To second order in perturbation theory, the new relationships for the Fermion operators 
are given by 



tl = «£°d£! + 4 6) C + 7^1 (4 a) cos k y sin((tf a - 9 b )/2) - cos k x sin 6> a )dg 



k+Q-a 



JhS v 

- (4 a) cos fc a sin # a + cos k y sin( (0 a - b ) /2)) 4+q -a | > ( 33 ) 

c£ = 4 b) 4 a 2 - 4 a) 4a + j^j -(4 a } cos ^ sin ^ + 4 b) cos K sin ((^ - ^)/2))4 a +q ,- Q 



+ 



(~4 a) cos k y sin((0 a - 6 )/2) + 4 b) cos k x sin^)4% _ I. (34) 



On the r sublattice, these equations imply that 9 r — ip r is related to the band and harmonic 
energies by 

n r 1 d E b 1 d E h 

- M9 r - = ___ = -j^^^ (35) 

which can also be obtained by minimizing the classical energy £ (9 a , 9b) = — JhS J2i n i cos(0j — 
ipi) + E h (9 a: 9 b ) with respect to 9 a or 9 b . Hence, the torque exerted by the electrons on the 
local moments opposes the tendency of the local moments to return to the angles that 
minimize E^. 

(r) (r)i 

In terms of the Fermion operators d^ and d^ that diagonalize H b , the full Hamiltonian 
can be expanded in a power series in 1/ v^S as H = Hq + Hi + H2 + ■ ■ ■■ The first-order term 
may be written 



v (r) v (s) a {b) 



where 



' k,q, r,s 

2t V ld {r)] d {s) (v {r) r (s) o(°) 4- v {r) v is) a^\ 

^/^y^Z^ j"k+q,T a k+Q,TV k+q X k «q + ^k+q^k «q ) 

+4 r |Ua4? (4 s) 4lq4 a) + W?) ) + h.c, (36) 

= cos ky sin((0 a — 9 b )/2) — cos k x sin 9 a , (37) 
= — 4^ cos k y sin((0 a — 9 b )/2) — cos k x sin 9 b , (38) 
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and x£ = — a^k+Q; Vk = vL+Qj V L = M k'i an d V L = ~ u \t > ■ The second term in Eq. (jHH|) 
is produced by the torque exerted by the electrons on the local moments. Notice that Hi 
is linear in the boson operators. The expectation value of the Fermion factor multiplying 
the boson operators in Eq. lfHKJ) vanishes provided that 9 r satisfy Eqs.(j2H} and (f2T)|) . The 
second-order Hamiltonian can be written as 

2 ~ N 2^ /^"^"k-qi.a^k-qa,/?] "qi "q 2 %-qi%-q 2 + "qi «q 2 %-qi U k-q 2 f" l JJ J 

k,qi,q 2 ,r,s a,/3 \ ) 

Since J^/JS and t/JS are of order 5°, Hi/JS 2 is of order 1/y/S and H 2 /JS 2 is of order 
l/S. 

To eliminate the first-order term in H and to express the Hamiltonian in terms of the 
rue SW operators for the total spin Sj )to t — Sj + Sj, we perform the unitary transformation 
H' = e~ u He u where U is constructed to satisfy [U,H ] = Hi. This transformation 
produces a modified second-order term H' 2 = H 2 + [U,Hi]/2. The anti-Hermitian operator 
U that fulfills these requirements is 

r/--2/m/^ V \d {r)] d {s) ( a ia \, {r) + A (r) v {s) \ - \ 

U - AJR\ 2^ ) "k+q,|"kT {% M k+q M k +«q %+q^k ) l 3 ) ( r ) ( 

" JV k,q,r,s I e k — C k+q — ZJ-^O ) 

It 

VSN k,q,r,s 

+ >)t A s ) ( n ( a hi {s) r (r) 7/ (r) \\ 1 he (AO) 

+ tt k+qa a k+Q,|V q U k+Q X k+q+Q + tt q "k+Qi/k+q+Q j L( s ) _( r ) ' iX -> l^ U J 



2* J j(r)t j(«) / (a) (r) (s) (6) (r) (s) 

2 < ] a k+q,t a k+Q,t V fl q M k+q X k ~+~ fl q V k+q^k 



C k+Q e k+q 



which is clearly of order {JS 2 )/\fS. 

Performing the unitary transformation and taking the limit J^S — > oo produces a rather 
complicated expression for H 2 . This expression can be considerably simplified by tracing 
over the Fermion degrees of freedom, thereby replacing the combination of Fermi operators 
^kt^k'T by ^ s expectation value 5 rs 5k,k'/(?k )■ Then H' 2 can be written as an effective 
Hamiltonian for the SW operators a^p and that takes the same form as the harmonic 
Hamiltonian of Eq. © but with revised coefficients A^'^ = A^'^ + (7 k + D^'^ and B k r ' s ^ = 
Bl r ' s) +4 ns) . The new contributions to these coefficients are 



Cl r ' s) = -jj £ /(e?) |< )2 (2f cos £ a cos(k + + e«)5 r ,A,a 
+vj) 2 (2t' cos6 b cos(k x + q x ) + eW)5 r , b 8 s>b 
+2t'cos((9 a - 6 b )/2) cos(k y + q v )u®v®(6 rta 5 atb + 5 r , b 5 s , a )\, (41) 
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(P>S ) _ At' JS 2 ^ / (I) (m) (J) (m) \ 

q,i,m 

C q+Q C k+q 

(r , s) _ 2t' 2 JS 2 v / (m) (J) (m) (0 \ 

-^k — V"q+Q X k+q+Q°r,a + ^q+Qi/k+q+Q^b^ 

q,i,m 

(«k+q4 + ^kiq^ >S s,b) ^ (i) -• (43) 

C q+Q C k+q 

Spin excitation frequencies are given by Eq. (jl3|) with A k and B^' s ^ replacing the 
harmonic coefficients. Note that the torque terms in Hi produce D^' s ^ and E^' s \ both of 
which vanish in the FM phase. In the CAF phase, these terms are required to preserve 
rotational symmetry and the relations ujq = and ujq = B. 

An analytic expression for the SW frequency is possible in the FM phase: 

co k = B + JS{3 - r) + {rj - 1) cos k x ) + K{2 - cos k x )/AS 

-JSy/(l + 7?) 2 (1 - cos k x ) 2 + (2 + K/AJS 2 ) 2 cos 2 k y , (44) 

which generalizes Eq. (115)1 for the Villain model. The FM phase becomes unstable when 
ujq = 0, which yields the same condition for the CAF-FM phase boundary as Eq.(|3U|). 
In the FM phase, the SW stiffnesses are the sum of DE and Heisenberg contributions: 
D x svi = K/8S + J3(l - rj)/2 and D'^ = K/8S + JS, both independent of field. For J = 0, 
these results agree with the SW frequencies of the DE model first obtained by Furukawa 
[ssT ]. When r) = —1, the SWs are isotropic because all nearest-neighbor interactions equal 
J. 

In the CAF phase, the SW frequency and stiffnesses are solved numerically by integrating 
over the first Brillouin zone to obtain C^' s \ -£? k , and E^' s \ Results for the SW frequency 
are plotted in Fig. 11 for rj = 3, B' = 0, and several values of t'. For all t' < t' c pa 21.2, the SW 
frequency vanishes at both the FM and AFM Bragg vectors. Above the phase separation 
region around t' pa 10.0 but below t' c , uj^ develops kinks that correspond to transitions across 
the neck of the a FS (Aqi = Q — ki pa 0.147rx for t' = 10.2) and across the length of the b 
FS (Aq 2 = Q — k 2 ~ 0.337TX for t' = 10.2). These transitions are sketched in Fig. 10. Since 
the nesting across Aq 2 is much weaker than across Aq 1; the kink in the SW frequency at 
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(0,0) 

(k x /7u,k y /7t) 




FIG. 11: The SW frequencies for p = 0.66, 77 = 3, B' = 0, and various values of t' . Different 
frequency scales are used on either side of k = 0. 

k2 is much weaker than the one at ki. As t' — > 10.0, Aq« — > so the kinks at ki and k2 
merge with Q. Whereas uj^ is a non-monotonic function of t' along the x direction, it is a 
monotonically increasing function of t' along the y direction. 

As plotted in Fig.6(b) for rj = 2, the SW stiffness in the x direction reaches a minimum 
at t' c ~ 13.2, above which both and are linearly increasing functions of t'. The 
stiffnesses in the x and y directions cross in the region of phase separation, where the SW's 
become isotropic. Notice that D™ increases by about 65% as t' increases from zero to t' c . 
For rj = 3, D™ doubles in this interval j^J 

With increasing rj, the SW stiffnesses plotted in Fig. 7(b) for t' = 10 and B = again 
cross in the region of phase separation around r] = 3.0. For larger 77, D* w > -D| w . Due to the 
strong increase in with 7], the average SW stiffness actually grows as the local moments 
become more non-collinear. Similar behavior was found when t = in Fig.2, where D™ was 
also found to be an increasing function of 77. 

The effect of a magnetic field is shown in Figs. 8(b) and 12 for rj = 2 and t' = 3. Absent 
from Fig. 8(b) is the region of very low fields, where the behavior of the SW stiffness is 
complicated by the dependence of the function D^ w (k a ,B) = (o; k — B)/k^ on the ratio 
29^ . For larger fields, there is a gradual increase in Z)^ w and a decrease in D^ w as 



K/VB' 

B grows from to B c w 4.SJS. The drop in D* w with field is caused by the growth of the 
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(k x /7t,k y /7t) 

FIG. 12: The SW frequencies for p = 0.66, rj = 2, t' = 3, and various values of B' . 

SW gap ujq = B at k = while the SW frequency continues to vanish at Q = (ir, 0) for 
B < B c , as seen in Fig. 12. Hence, the average SW stiffness drops from roughly 1.2JS at 
B = to 0.54JS at B c . Because a magnetic field very quickly destroys the region of phase 
separation, the crossing point of the x and y stiffnesses does not coincide with a region of 
phase separation. As expected, the SW stiffnesses are independent of field above B c . Notice 
from Figs. 11 and 12 that the SW velocity at Q softens within the CAF phase as t approaches 
t c or B approaches B c . This behavior is quite general and occurs at the CAF-FM phase 
boundary of any magnetic system jjnj]. 



V. DISCUSSION AND CONCLUSION 



Applying the results of this work to the manganites requires an estimation of the relevant 
parameters. Using S = 3/2, r/JS ~ 8.6 meV, i] ~ 1 [jifj], and t « 200 meV for the metallic 
phase gives t' ~ 32 and B' « 0.0285(T). When rj = 1, the critical values of the DEV model 
are t' c = 5.3 in zero field and B' c = 2.5 for t' = 0. Hence, the metallic value for the hopping 
energy is more than sufficient to stabilize the FM phase but an applied field of 3 T cannot 
align the local moments when t' is much below t' c . 

Metal-insulator transitions in many of the man gan ites have been successfully interpreted 
within the framework of percolation theory Issj]. When the FM fraction / exceeds 
a critical value f c , a FM backbone spans the sample and the system becomes metallic. 
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The percolation threshold f c in those manganites lies between 0.17 and 0.19. But it is 



rather difficult to explain the jump in the SW stiffness observed in Pro.7Cao.3Mn03 [11 1 
within the framework of percolation theory. Assuming that the metal-insulator transition 
in Pri_2.Ca2.MnO3 is produced by the percolation of metallic FM clusters [10], there is no 
reason why the SW stiffness should change when a FM backbone spans the sample. Within 
percolation theory, one would more naturally expect the SW stiffness to gradually increase 
with field as a growing fraction of the sample becomes metallic. 

Estimates of the FM fraction / in Pr 1 _ 2 ,Ca ;c Mn03 manganites also pose a difficulty for 
percolation theory. The accepted range for the percolation threshold f c within continuum 
models falls between 0.15 and 0.22. While one experimental estimate of f c ~ 0.08 at the 
metal-insulator transition lies well below this accepted range, a recent estimate of f c ~ 0.6 
Il3l based on more complete magnetization data lies well above this range. At zero field 

n q 

and low temperatures, estimates for / range from 0.3 [13J to 0.5 |9J]. If those estimates are 
correct and the FM regions are metallic, then the percolation threshold would already be 
exceeded and the low-temperature phase would be metallic in zero field. 

On the other hand, our results clearly indicate that the jump in the SW stiffness at a field 
of 3 T cannot be produced by simply aligning the AFM regions while keeping the bandwidth 
~ t fixed. As found in Fig.8(b), a magnetic field would act to suppress rather than enhance 
the average SW stiffness under those conditions. Moreover, a field of 3 T corresponding to 
B' 0.086 <C B' c will have little affect on the alignment of the local moments. Whereas 
experiments indicate that the AFM regions begin to shrink under a magnetic field greater 
than 3 T, a sizeable fraction of the sample remains AFM even above the metal-insulator 



transition 
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13|. 



The large increase in the SW stiffness predicted by our work as if increases from to 

t' c suggests an entirely different scenario: at the critical field, electrons in the insulating 

and possibly canted FM regions delocalize as the hopping energy t sharply increases. Since 

the integrated optical weight is just K ~ t, the jump in the hopping energy at 3 T should 

rn 

be observable in optical measurements. Indeed, measurements by Okimoto et al. |7j on 
Pr .6Cao.4Mn03 reveal a rapid rise in (t(u) and a rapid drop in the CDW gap near the 
critical field. If the percolation threshold for the FM regions is exceeded when the electrons 
delocalize, as suggested by the measurements discussed above, then the jump in the SW 
stiffness will coincide with the metal-insulator transition. Otherwise, the metal-insulator 
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transition will occur at a slightly larger field. 

Supporting this picture are the plethora of probes that can produce a metal-insulator 
transition in Pr 1 _ a; Ca x Mn03 (0.3 < x < 0.4). Besides a magnetic field, application of an 



electric field j^, high pressure j40||, exposure to x-rays 4l|, and exposure to visible light 42 1 



all induce a metal-insulator transition. The identical resistivities produced by application 



of a magnetic field or exposure to x-rays 



41] suggest a common mechanism: the excitation 



of charge carriers out of polaronic 
subsequent relaxation of the lattice 



;raps produced by the electron-lattice coupling. The 
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4l| may prevent retrapping of the electrons. 



, which gives the lifetime r(k) = 
43j of the spin-diffusion coefficient 



The jump in the SW stiffness as a function of field in Pro.7Cao.3Mn03 also bares 
a striking resemblance to the jump in the SW stiffness observed in Lai_2.Caa.MnO3 at the 
metal-insulator transition with x ~ 0.22 4sj|. In fact, the sizes of the SW stiffnesses on either 
side of the doping-induced transition in Lai-zCa^MnOs are almost exactly the same as on 
either side of the field- induced transition in Pr .7Cao.3Mn0 3 . This suggests that there are 
also large increases in the hopping energy and electronic kinetic energy in Lai-^Ca^MnOs 
at the critical concentration of x ~ 0.22. 

Compare the jump in the SW stiffness across the metal-insulator transition with the 
smooth dependence of the spin-diffusion coefficient A 4], [ 
of paramagnetic spin fluctuations. Measurements 
were performed just below Tq, where paramagnetic spin relaxation occurs within polaronic 
regions of the FM phase. While _D SW is a linear function of the electronic kinetic energy in 
the FM phase, Ax depends only on doping and is independent of the electronic bandwidth 
W in the low-temperature limit T <C W [4]. So long as the bandwidth of the polaronic 
regions remains large compared to the temperature, the spin- diffusion coefficient will not 
change across the metal-insulator transition. 

Except in a very narrow range of parameters, phase separation is absent in the DEV 
model. Phase separation appears quite commonly in DE models with AFM Heisenberg in- 
teractions and even occurs near p = 1 in a DE model without Heisenberg interactions 
I. Results for the DEV model contrast with both a DE model with AFM interactions 



id a 

44]. 



between all neighboring local moments 
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22. 



and a DE model with FM interactions 
in plane but AFM interactions between neighboring planes 2] • Due to the high symmetry 
of both models, AFM order is not frustrated when t = 0. But in both cases, phase separa- 
tion occurs before the AFM interactions become strong enough to cant the spins. Work by 
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Golosov 



generacies 



231 ] indicates that the canted phase is destabilized by the presence of local de- 



24j that are absent in the DEV model. Our results suggest that phase separation 



only occurs within a very narrow range of parameters when the AFM Heisenberg interac- 
tions are frustrated and local degeneracies are absent. For doping concentrations away from 
multiples of a quarter filling, long-range orbital ordering is impossible and magnetic frus- 
tration may be present in a wide range of manganites. Therefore, the DE model may not 
provide as straightforward a pathway to phase separation in the manganites as believed. A 
much more important role may be played by the quenched disorder associated with chemical 
inhomogeneities 

Roman and Soto j^| pointed out that the nature of the FM and AFM regions can be 
probed by measuring the SW spectrum about the AFM Bragg vector Q. In a CAF, the 
SW's are not gapped for any field B below B c , as shown in Figs. 4 and 12. By contrast, the 
SW's of a collinear AFM contain a gap at Q that grows linearly with field. So measurements 
of the SW spectrum about Q may be used to determine whether local regions are canted. 
As discussed elsewhere [29] and seen in the inset to Fig.3(b), the rapid increase in the SW 
stiffness for fields near B* is another signature of a canted phase. The absence of this 
behavior in Prx.a-Ca^MnOa for 0.3 < x < 0.4 strongly suggests that the local AFM regions 
are not canted. 

Since CE-type AFM ordering occurs at half-filling, the local charge ordering observed in 
the manganites Pro.6Cao.4Mn03 and Lao.7Cao.3Mn03 would be simplified if the 

polaronic regions were rich in holes and poor in electrons compared to the bulk. Our model 
provides a natural explanation for this behavior, since the electronic fraction on b sites is 
substantially smaller than the fraction on a sites as the electrons avoid regions with more 
pronounced AFM order. 

This paper has studied the general effects of AFM interactions and non-collinearity on 
the magnetic ordering and spin dynamics of DE systems. The competition between DE 
and AFM interactions is responsible for several interesting properties. In a narrow region 
of hopping energies, weak phase separation occurs as the FS topology changes from closed 
to open. Because electrons prefer to occupy sites that are coupled by FM interactions, a 
CDW appears in the absence of a CDW gap. For finite Hund's coupling, the electron spins 
are more closely aligned than the local moments of the CAF. The CAF becomes unstable 
above critical values of field and hopping energy and below a critical value of rj. Perhaps 
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the most surprising result is that DE changes none of the qualitative features of the CAF 
state. The SW gap u at k = remains equal to B and the gap ujq at the AFM Bragg 
vector Q continues to vanish for any hopping energy below t c . The average SW stiffness 
still softens with magnetic field and hardens with strengthening AFM interactions; the SW 
velocity at Q still vanishes at the CAF-FM transition. Clearly, a great deal can be learned 
about the general properties of CAF's, even of the more complex variations that appear in 
many manganites, by considering simple models such as the DEV model. 
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